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Abstract 

Present and future expected limits on interactions between dark matter 
and various quarks are thoroughly investigated in a model-independent way. 
In particular, the constraints on the interactions from the Large Hadron Col- 
lider (LHC) experiment are carefully considered with focusing on mono jet + 
missing transverse energy {jPt)i mono b-jet + p?>, and top quark(s) + $x 
channels. Model-independent upper limits (expected limits) on the cross sec- 
tion times acceptance for non-standard model events are derived at 7 TeV (8 or 
14 TeV) running of the LHC experiment. With assuming that the dark matter 
is a singlet real scalar or a singlet Majorana fermion, we also put constraints 
on several operators describing its interactions with up, down, strange, charm, 
bottom and top quarks. These constraints are compared to those obtained by 
cosmological and astrophysical observations of the dark matter. 



1 Introduction 



The origin of the dark matter in the Universe and its nature are the most outstanding 
mysteries in modern particle physics, cosmology, and astrophysics. Though the prob- 
lem is not resolved yet, many theoretical and experimental efforts have been made to 
search for the dark matter, and those give precious information about several dark 
matter interactions. Theoretically, the weakly interacting massive particle (WIMP) 
is one of the promising candidates for dark matter [1] and the couplings of the in- 
teractions involving the WIMP are predicted to be of the order of the electroweak 
(EW) couplings. Experimentally, direct detection experiments of dark matter put 
very severe limits on the scattering cross section between the dark matter and a nu- 
cleon [2], while indirect detection experiments give constraints on self-annihilation 
cross sections of the dark matter into various final states [3]. 

Among dark matter interactions, those to quarks and gluons are particularly well 
investigated thanks to the Large Hadron Collider (LHC) experiment. Importantly, 
the LHC experiment enables us to explore not only the interactions relevant to 
direct and indirect dark matter detections but also the interactions irrelevant to 
those detections. According to this fact, many studies have already been performed 
assuming a definite relation among the dark matter interactions with various quarks 
or using an analysis with simple collider simulations [I]. In this article, we thoroughly 
investigate present and (near) future expected limits on the dark matter interaction 
with each quark with focusing on mono jet + missing transverse energy {$t)i mono 
b-jet + p T , and top quark(s) + fi T channels, which are all derived from dedicated 
collider simulations with use of a model-independent method. 

In the next section (section [2]), we first summarize dark matter interactions as- 
suming that the dark matter particle is a real scalar or a Majorana fermion which is 
singlet under SM gauge groups. These dark matter candidates are frequently used 
in many literatures and can be regarded as the simplest example of a dark matter 
particle. We then put limits on the operators from cosmological and astrophysical 
observations of the dark matter in section [3J We next consider what kinds of dark 
matter signals are generally expected at the LHC experiment in section H] and give 
model-independent upper limits (expected limits) on the cross section times accep- 
tance for dark matter signals at 7 TeV (8 or 14 TeV) run. It turns out that these 
limits are very useful to constrain any dark matter quark interactions. Using the 
obtained limits, we also explain how to put a limit on each operator (discussed in 
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section |2]) describing a dark matter interaction with various quarks. We finally con- 
sider how severely the LHC experiment can put limits on those operators in section 
and compare the limits with those obtained from the cosmological and astrophysical 
observations. Section [6] is devoted to a summary of our discussions. 



2 Dark matter interactions 

Dark matter interactions with various quarks are discussed in this section assuming 
that the dark matter is a real scalar or a Majorana fermion which is singlet under 
SM gauge groups. We introduce operators describing interactions between the dark 
matter and SM particles up to the mass- dimension of six. Decoupling and/or weakly 
coupled heavy physics behind the interactions are implicitly assumed, so that the 
whole theory is renormalizable. Because of the stability of the dark matter, which 
is guaranteed by imposing a Z2 symmetry, the operators have to involve two dark 
matter fields. The real scalar dark matter is denoted by <fr, while the Majorana 
fermion dark matter is denoted by x hi following discussions. 

The effective lagrangian for the scalar dark matter field <p is simply given by 
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where £sm is the SM lagrangian and is the invariant mass term of the dark 
matter. Since several interactions will contribute to the dark matter mass after the 
electroweak symmetry breaking, we use the notation to describe the physical 
dark matter mass throughout this paper. The mass dimension of operator is 
denoted by n, while A is the cutoff scale. The effective lagrangian describes the 
physics below the scale. The complete set of higher dimensional operators, which 
are invariant under both Lorentz group and SM gauge groups, is given by 
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where H, L = (v, ei) T , E, Q = (ul, ^l) t , D, and U are the higgs doublet, the lepton 
doublet, the charged lepton singlet, the doublet quark, the down-type singlet quark, 
and the up-type singlet quark, respectively. On the other hand, (V^„ = e^ va pV a ^) 
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is the field strength tensor of a gauge boson, namely, V^ v = B^, W^ u , and G" 
are those of the hyper-charge gauge boson, the weak SXJ(2) L gauge boson, and the 
gluon, respectively. The covariant derivative is denoted by D^. Operators 0\£ HE , 
O^Q HD , and O^Q HU are in fact flavor- dependent such as <p 2 QiHDj. In order to avoid 
dangerous flavor changing processes, we only consider flavor diagonal parts of the 
operators. This means that we are implicitly assuming that heavy physics has some 
mechanism to suppress flavor off-diagonal parts of the operators. 

Since we are interested in interactions between dark matter and various quarks, 
we are focusing on O^Q, HD . and O^Q. HU . (i = 1,2,3) among the operators listed 
above. After the electroweak symmetry breaking ((0\H\0) = (0,v) T /\/2 with v 
being about 246 GeV), these operators give following effective interactions, 

2 ^ 

A ( nt = 5Z [ Ui ( C 0oUtA + ^LhuJs) U i + di [c^q.hd. + ^Q^D^ di , (3) 

where cf Q . HU . {cfy. Hu .) and c$. HD . (c% iHD ) are real (imaginary) parts of the 
coefficients c (j>QiHDi and c^ Q . HUi in front of the operators 0^. HD . and 0^. HU .. 

As in the case of the real scalar dark matter field 0, the effective lagrangian for 
the Majorana fermion dark matter field x is a l so given in the same form, 
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where we adopt the four component notation of the Majorana field x = X c with the 
superscript V being the charge conjugation. Invariant mass term of the x field is 
denoted by M x , which is taken to be real by an appropriate redefinition of x- We use 
the notation m x to describe the physical mass of the dark matter x- The complete 
set of effective operators involving two dark matter fields is given by 

0% = (XX)\H\\ O x % = (xri.xWl.U) 

0% = *(X75X)|# I 2 , Of D = {X1^X){D1,D) 

Of H = (xri 5 x)(HU^H), 0f L = ( X r^x)(L^L) 

o { 2 = (xri 5 x)(Qi,Q), of E = (xYkxXe^e). 

Because of the same reason as the scalar dark matter case, we focus on flavor diagonal 
parts of the operators O^q, O^, and 0^% which give following interactions, 

Ant = tt~ ( c xu. p n + c xQ, P l) Ui + da„ (c xDi Pr + c xQi P L ) di] , (6) 



i=i 



where all coefficients c x u v c x d and c X Q i are real numbers. 
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3 Cosmological and astrophysical constraints 



We are ready to put limits on each operator describing interactions between the 
dark matter and various quarks. In this section, we consider various limits obtained 
from cosmological and astrophysical observations. The current and future expected 
limits on the operators obtained from the LHC experiment will be discussed in 
the next section. For the real scalar dark matter 0, following twelve operators, 
4> 2 UiUi, (f> 2 Uiij5Ui, (p 2 didi, and 4> 2 dii^^di are considered, where index % runs from one 
to three. On the other hand, for the Majorana fermion dark matter Xi following 
nine operators, (xYl^x)^^^), (x7 m 75X)(^7m^r^), and (xYl5X)(u i >y fl P L u i + 
dci^Phdi) are considered. Because new physics model generally involves several 
operators simultaneously with appropriate coupling constants, putting limits on each 
operator should be regarded as a toy model. With keep in mind this fact, we simply 
assume that the effective lagrangian has only one of the operators mentioned above 
and explain how to put a limit on each operator in this and next sections. 



3.1 Cosmological limit 

A cosmological limit is obtained from the thermal relic abundance of the dark mat- 
ter [5J. When the annihilation cross section of the dark matter is smaller, the resul- 
tant abundance is larger, leading to the over-closure of the universe. We therefore 
obtain a lower limit on the coefficient of the operator from the cross section. 

For the scalar dark matter <p having interactions (v(p 2 / \/2A 2 )(c^qq + c^qij^q), 
where q represents a quark (q — v,i, di) with m q being its mass, the annihilation cross 
section of the dark matter <fi (times relative velocity v Te \) is evaluated to be 

3v 2 
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(c^y (1 - m)lm\fl 2 + (c^Y (1 - mj/mj) 1 ^ + Q(e 



where e = (s — 4m|)/(4m|). Using a semi-analytical formula for the dark matter 
abundance [6], we obtain a lower limit on the coefficient vc^) /A or \J /A. 

On the other hand, for the Majorana fermion dark matter x having a interaction 
(2c/A 2 )(x7 M 7 5 x)(g7 A1 P R ( i )g), the annihilation cross section is given by 

av Iel \ xx ~ i^L (4m 2 e + 3mJ/2) yjl - m 2 /m 2 x , (7) 

where e = (s — 4m 2 )/(4m 2 ) again. As in the case of the scalar dark matter, we 
obtain a lower limit on the coefficient yfcj A using the same semi-analytical formula. 
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3.2 Astrophysical limit 1 (Direct detection) 

Direct detection experiments of dark matter put a stringent limit on each opera- 
tor when the operator contributes to the spin-independent scattering cross section 
between the dark matter and a nucleon. In this paper, we adopt the result of the 
XENON100 [2]. For the scalar dark matter 0, as shown in equation (131) , there are two 
types of interactions; scalar and pseudo-scalar interactions. Since the pseudo-scalar 
interaction is not constrained, the operators of the scalar dark matter are limited 
only through the scalar interaction. 

Using the formula discussed in reference [7] , the spin-independent scattering cross 
section between the dark matter and a nucleon is given by the following formula when 
the dark matter couples to light quarks (w, d, and s quarks); 

0"si " 



27rA 4 (m 9 j + m^m 1 



where — 940 MeV is the nucleon mass, while fx q is the parameter determined 
by the hadron matrix element (N\qq\N). The values of the parameters for u, d, and 
s quarks are f Tu — 0.028, fad — 0.028, and f Ts ~ 0.009, respectively [8]. On the 
other hand, when the dark matter couples to heavy quarks (c, b, and t quarks), the 
formula of the spin-independent scattering cross section is changed to be 

{c^) 2 2v*m%P TG 
aSI 729nA^m (j} + m N ) 2 m 2 q ' [ ' 

The parameter Jtg is determined by the hadron matrix element (N\G a ^ v G a ^ v \N) 
and its value is evaluated to be ]tg — 0.9431 through the trace anomaly relation, 
namely fr u + frd + fr s + frc = 1- The essential reason why the scattering cross 
section depends on the gluon hadron matrix element is that the dark matter can 
be scattered with a nucleon through the interaction with gluons inside the nucleon, 
whose process is mediated by the one-loop diagram of the heavy quark. 

On the other hand, for the Majorana fermion dark matter x, all interactions with 
quarks come from the axial-vector interaction of the dark matter, namely X1^15X- 
Only the spatial components of the interaction are survived at the non-relativistic 
limit, so that the corresponding quark interaction becomes qYl5Q- These operators 
therefore do not contribute to the spin-independent scattering cross section but con- 
tribute to the spin-dependent one. Since a limit from direct detection experiments 
through the spin-dependent interaction turns out to be much weaker than that from 
the LHC experiment, we do not consider the constraint in this article. 
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3.3 Astrophysical limit 2 (Indirect detection) 

Another astrophysical limit is obtained by indirect detection experiments of dark 
matter. Several products from dark matter annihilations at the galactic halo, for 
example, gamma-rays, positrons, electrons, anti-protons, and anti-deuterons, are 
utilized to detect the dark matter. Among those, one of the reliable limits on the 
annihilation cross section of dark matter comes from the indirect detection exper- 
iment through gamma-ray observations. In particular, the observations from the 
galactic center and milky- way satellites currently give the most stringent limit. In 
this article, we put a limit on each operator according to reference [9], in which 
a limit on the dark matter annihilation cross section obtained by the gamma-ray 
observation from the galactic center at the Fermi-LAT experiment is presented. 

Differential gamma-ray flux in a direction ip originated in dark matter annihila- 
tions at the region around the galactic center is given by the following formula, 



where a is the annihilation cross section with v m \ being the relative velocity between 
incident two dark matters, while tbdm is the dark matter mass. The mass density of 
dark matter (the dark matter profile) at the position I is denoted by p(l). The square 
of the dark matter density profile is integrated over along the observed line-of-sight. 
In our analysis, we have used the "Cored (R c = lkpc)" dark matter profile, which is 
one of the profiles used in the reference and gives the most conservative limit on the 
annihilation cross section. The fragmentation function dN^/dE, which is nothing 
but the gamma-ray spectrum produced per annihilation, has been estimated using 
the PYTHIA code [IT] for various final states like uu, ss, bv, ti, etc. 

Using the observational result of the Fermi-LAT experiment, model independent 
upper limits on the quantity (av vc \/m^ M ) x J^f dE[dN^/dE] for various energy bins, 
(Ei,E f ) = (0.3, 1), (1, 3), (3, 10), and (10, 100) in GeV unit, are presented in 
the reference. The upper limit on the annihilation cross section, av ie \, is therefore 
obtained when the dark matter mass and the fragmentation function are given. 
Because the dark matter velocity at the present universe is three orders of magnitude 
smaller than the speed of light, the annihilation cross section is expressed by the s- 
wave component only. For the scalar dark matter, this is expressed by 
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where the interaction of the dark matter is assumed to be (v 2 /v / 2A 2 )(c^gg + 
c^qi^q). On the other hand, for the Majorana fermion dark matter x, the s-wave 
component of the annihilation cross section is given by the following formula, 

6c 2 m 2 f / 

The dark matter is assumed to have the interaction {^c/ A 2 )(xifil5X)(Ql fl PR(L)q)- 
Since the annihilation cross section is proportional to the fermion mass squared be- 
cause of the helicity suppression, only the interactions involving the second and third 
generation quarks are constrained by the indirect detection experiment. 

In this section, we have discussed the method to put limits on the dark matter 
operators using cosmological and astrophysical observations. Obtained limits on each 
operator are summarized in figure [3] for the real scalar dark matter and figure H] for 
the Majorana fermion dark matter in section where limits from cosmology, direct 
detection, and indirect detection experiments are shown as brown, green, and blue 
lines in both figures. Physical interpretation of the results are discussed in section [5] 
in comparison with the limits obtained from the LHC experiment. 

4 Constraints from the LHC experiment 

We first discuss possible signals of dark matter at the LHC experiment. Dark mat- 
ter particles are pair-produced from quark-quark interactions through higher dimen- 
sional operators discussed in section [2j However, the mere pair production does not 
produce missing transverse energy an d cannot be seen at the experiment. For 
this purpose, we need to pair-produce dark matter particles in association with jets 
or photons. In this case, the dark matter pair recoils against the visible particles 
and produces $ T . The generic signature of dark matter production at the LHC 
experiment is thus $ T plus a few number of jets or photons. 

In view of collider searches, we can divide the dark matter effective operators 
into three parts, i.e., the operators containing light quarks, bottom quarks and top 
quarks respectively. In the first case, we can pair-produce dark matter particles with 
gluon/quark radiations from initial quark/gluon legs and thus observe the signal in 
the mono-jet plus $ T channel. In the second case, we need two bottom quarks in the 
initial state that may come from the gluon splitting. This means, we may produce 
dark matter pair plus b jets. The dark matter interaction with top quarks is very 
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Figure 1: $j> distributions for the two operators described in the text at 8 TeV run. 

different from two other cases. Here we need to produce two dark matter particles in 
association with two top quarks and the signal looks like top pair plus $ T . This final 
state is particularly interesting as it comes from different beyond the SM processes. 

From above discussions, it is clear that the signals will be visible as events with 
P T plus a few number of SM particles and the shape of the fJ T distribution may 
be used to separate out those from SM backgrounds. However, unlike the R-parity 
conserving SUSY, $ T comes from the recoil of jets against dark matter particles and 
thus produces featureless rapidly decreasing distribution for dark matter coupled to 
light and bottom quarks. In case of top quark and dark matter effective operators, 
situation is slightly different as $ T can also come from semi-leptonic decays of W 
bosons. It is therefore important to compare $ T distributions for these two cases. 
The $ T distribution becomes harder with increasing the dark matter mass although 
the overall production cross section decreases with the mass. 

For illustrative purpose, we take the two operators (2x7 m 7 5 x)(m7 a1 P r m)/A 2 and 
(2x7^75%) (t-f^Pnt) / A 2 assuming m x = 100 GeV and A = 1 TeV. The fJ T distribution 
is shown in figure [T]for 8 TeV run with 15 fb _1 integrated luminosity. It can be seen 
from the figure that, in case of top quark dark matter effective interaction, the $ T 
distribution is slightly harder compared to the light quark case although the cross 
section is extremely small. For this reason, it is not possible to probe such dark 
matter that couples very weakly to top quarks from current 8 TeV run and we have 
to consider high luminosity LHC run (100 fb" 1 ) with a/s = 14 TeV. 
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4.1 Model independent upper limits 

We now discuss model-independent upper limits on the cross section times accep- 
tances for new physics interaction involving dark matter with three channels: mono 
jet + $ T , mono b-jet + $ T and top quark(s) + $ T searches. In order to estimate 
the prospects with 8 and 14 TeV runs of the experiment, we perform a Monte-Carlo 
(MC) simulation. Madgraph5 [10] is utilized to generate parton level events and the 
events are interfaced to PYTHIA6.420 [11] to deal with showering and hadroniza- 
tion. For detector simulation, we utilize Delphes2.0.2 [12] and resolution parameters 
are chosen based on the CMS detector performance given in reference [T3] . 

4.1.1 Mono jet + $ T 

The upper limit on the non-standard model contribution to the mono jet + $ T 
channel have already been reported in references [T3"j dU [TJ] for 7 TeV run. In this 
paper, we take the results of reference [H] and set the constraints. We also estimate 
the prospect for 8 TeV run by the MC simulation imitating the search discussed in 
the reference. The basic selections used in the analysis are the following: 

1. Event has the missing energy 

t T > 200 GeV. 

2. The leading jet has pt >110 GeV and \r)\ < 2.4. 

3. The number of jets with p? > 30 GeV is smaller than 3. 

4. There is no isolated leptons or tracks with p^ > 10 GeV. 

5. The second jet has pt < 30 GeV or the difference of the azimuthal 
angle between the leading and the second jets is smaller than 2.5. 

Here, notations pt and 77 denote the transverse momentum and pseudo-rapidity 
respectively. In reference [T4], with additional cuts $ T > 250, 300, 350, or 400 GeV 
applied, the upper limit on the number of the non-standard model events has been 
obtained. We express the results in terms of the cross section times the acceptance 
(a X A), which are summarized in the table shown below; 



Cut on $ T (GeV) 


250 


300 


350 


400 


a x A (fb) 


120 


73.6 


31.6 


19.0 



In order to estimate the prospect with the 8 TeV run as accurately as possible, 
we first perform the MC simulation with 7 TeV run and compare the background 
numbers with those of reference [H] and obtain normalization factors required to 
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$ T > 350 GeV 


2065 


4674 


28 


$ T > 400 GeV 


928 


2337 


20 


$ T > 450 GeV 


473 


1275 
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$ T > 500 GeV 


223 
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Table 1: Cut How table for the SM backgrounds at the mono jet + $t channel. The 
center of mass energy of 8 TeV and the integrated luminosity of 15 fb^ 1 are assumed. 

adjust the difference. The same normalization factors are used to calculate the 
background for 8 TeV run. Note that the cross section obtained by MadGraph and 
PYTHIA is only at LO, and one needs K factor to obtain the correct normalization. 
Instead using the K factor, we compare our simulation to the distribution in the 
reference. As is shown in the reference, W — >■ Iv + jets, Z — > vv + jets and ti + jets 
are the dominant backgrounds, where I means the charged leptons. Other possible 
sources, QCD, single top, etc. are negligible. The cut flow for the standard model 
events is shown in table [1] assuming the integrated luminosity of 15 fb _1 . 

We have calculated the expected 95% C.L. upper limit on the non-standard model 
cross section times the acceptance, which is summarized in the table below. Here, 
we have assumed that the systematic uncertainty associated with the estimation of 
the SM background events is 10%, which is same as the one in reference [T4] . 



Cut on $ T (GeV) 


250 


300 


350 


400 


450 


500 


a x A (fb) 


459 


191 


89.1 


43.6 


23.7 


13.7 



4.1.2 Mono b-jet + $ T 

The mono b-jet + $ T channel is expected to be useful to search for the dark matter 
which mainly couples to bottom quarks. In addition to the selection cuts used in 
the previous subsection (the cuts 1-5), we also require the following one: 

6. The leading jet is b-tagged. 
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W -> lu+jets 


Z — > vv+iets 


tt+jets 


$ T > 350 GeV 


33 


111 


23 



Table 2: Expected number of the SM events in the mono b-jet + $x channel. The center 
of mass energy of 8 TeV and the integrated luminosity of 15 fb" 1 are assumed. 

For the b-tagging method, we assume the following properties: 

• The b-tagging efficiency is 0.6 for real b jets. 

• The miss-tagging rate for light jet is 0.004. 

• The miss-tagging rate for c-jet is 0.1. 

These tagging efficiencies can be achieved in the CSVM tagger calibrated in refer- 
ence [IE]. In the case of the mono b-jet + fJ T channel, we should not raise the cut on 
$ T as strong as that of the mono jet + $ T channel. For large p T events, we expect 
large pt for b-jets, which lower the efficiency of the b-tagging. For this reason, we 
just fix the cut on $ T to a moderate value (>350 GeV) and set the constraint. 

The estimation of the number of the SM events after the cuts 1-6 and p T > 
350 GeV is shown in table [2] for 8 TeV run assuming 15 fb _1 data. Here, "jets" 
also include the jets originate in heavy flavors. Since the mono b-jet + $ T channel 
is almost similar to the mono jet + $ T channel, the dominant backgrounds in the 
previous subsection are also the dominant ones here. We use the same normalization 
factor as the one in the mono jet + $ T channel to estimate the number of the events. 

The expected 95% C.L. upper limit on the cross section times the acceptance 
for the non-SM event is estimated to be 2.9 fb. Here, we have again assumed that 
the systematic uncertainty in the estimation of the SM events is 10%. One may 
expect that the use of tighter b-tagging condition efficiently reduces the W and Z 
backgrounds. However, it turns out that half of the remaining W and Z events 
include the real b-jets with the b-tagging properties we adopt. Therefore, tightening 
the b-tagging condition does not improve the result of our analysis. 

4.1.3 Top quark(s) + $ T 

Among the SM quarks, top quark needs separate treatment. It is heavy and not 
contained in protons. If the dark matter couples mainly to top quarks, dark matter 
production cross section is highly suppressed in comparison with the case with light 
quarks. In this case, the dark matter must be pair-produced associated with two top 
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quarks, and a top quark decays to hard b-jet, lepton, and neutrino that contributes 
to additional missing momentum. Therefore a specialized search strategy is required 
for such scenario. In this paper, we concentrate on the prospect for 14 TeV run 
because the production cross section is too small at 7 and 8 TeV runs. 

The dominant background process to the top quarks + $ T channel is the ti 
production. The difference between the background process and the dark matter 
production in association with top quarks is that there is another source of fJ T in 
addition to neutrinos from the leptonic decay of W boson. In order to suppress the 
SM ti events, we consider zero lepton mode in which the $ T from leptonic decay of 
W boson is suppressed. We therefore impose the following basic selection cuts: 

1. There is neither isolated electron nor muon with p? > 20 GeV. 

2. There are at least four jets with p T > 100, 80, 50, 50 GeV 
and one of them has to be b-tagged. 

3. The azimuthal angles between the missing transverse 
momentum and the leading four jets are larger than 0.2. 

4. M cff > 1000 GeV and $ T /M cE > 0.3. 

M e ff is the scalar sum of the transverse momenta of the leading 4 jets and $ T . 

Since the top quark(s) + $ T channel resembles the supersymmetry search with 
multi-jet + $ T [17], the dominant backgrounds are again W —> Iv + jets, Z — > vv + 
jets and ti + jets. Interestingly, we find that My variable [18] is useful to separate 
out signals from SM tt events. In order to define the variable, we first define two mega 
jets. If n jets are present in the event, there are 2"" 1 possibility for reconstructing 
two mega jets. Among them, we choose the one which minimize the sum of the 
square of the two mega jet masses. From the four momentum of the two jets p? 1 , pP 2 
and the transverse missing momentum fJ T , is defined by 



For the SM ti events, it is expected that distribution has a cut off around m t , 
while that of the signal is expected to be broader. We show M|? distributions for 
standard model ti and for dark matter effective operator (2x7^75%) (ij^Pjit) /A 2 in 
figure H] assuming m x = 100 GeV and A = 1 TeV after imposing the selection cuts 
1-3. Based on this distribution, we impose additional cuts > 300, 400, 500, 
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Figure 2: Normalized distribution of the variable for the SM background (ti) and for 
dark matter effective operator (2%7^75%) (i-TuPRi) I 'A 2 events after imposing the selection 
cut 1-3. Center of mass energy of the LHC experiment is set to be 14 TeV. 



600, 700, 800 and 900 GeV. The cut flow for the SM backgrounds, namely ti +jets, 
W — > Iv + jets, and Z — > vv + jets backgrounds are shown in table [3j 

As expected, the ti+ jets background reduces drastically with cuts applied, 
while Z and W backgrounds are not reduced significantly by this cut. Here, the 
number of the ti events is normalized assuming the total cross section of 920 pb [19] . 
The numbers of the W — > Iv + jets and Z — \ vv + jets events are normalized 
with the matched cross section obtained from the MadGraph-PYTHIA package. We 
have calculated the expected 95% C.L. upper limit on the non-standard model cross 
section times the acceptance, which is shown in the table below. Here, the systematic 
uncertainty in the estimation of the standard model events was assumed to be 10%, 
which is typical of the supersymmetry search with hadronic mode. 



Cut on Mj: (GeV) 


300 


400 


500 


600 700 


800 


900 


a x A (fb) 


12.2 


10.2 


4.57 


1.87 0.839 


0.418 


0.208 



In this paper, only the result for the hadronic mode has been shown. On the other 
hand, we also have estimated the prospects using the one (leptonic) and two lepton 
(semi-leptonic) modes. For the leptonic mode, we use the My 2 variable constructed 
from $ T and the transverse momentum of two leptons. For the semi-leptonic mode, 



13 





tn^jeis 


/ V. Ill 1 1 T nf o 

Z; — ? t/t/-rJcLo 


t tt/-|-jcio 


fl„f 1 A 

OUT 1-4 


( 40 




1 IP. 
1(0 


iKZy -> OUU Vorti V 


^79*3 
( Zo 


9^ 
zoo 


1 7fi 
1 ( 


/i/T^ \ /inn p D \/ 
.My -> 4uu ^orev 


4 I Oo 




1 7/I 
1(4 


Ivlrp ^> OUU VJTC V 


1 Q^4 


9nn 

zuu 


1 41 

±4: J. 


My > 600 GeV 


659 


143 


95 
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50 
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94 


47 


19 


My* > 900 GeV 


21 


21 


14 



Table 3: Expected cut flow for the SM backgrounds in the top + Pt channel. The center 
of mass energy of 14 TeV and the integrated luminosity of 100 fb" 1 are assumed. 

the My variable constructed from $ T and the transverse momentum of a lepton 
are used to reduce the background. These methods are discussed in reference [20] , 
and those would be useful for the scalar top search. We find, however, that the 
constraints on the effective operators involving top quarks, which will be discussed 
in the next section, are weaker than the one with the hadronic mode. 

We are now ready to discuss how we can put the current and future expected 
limits on each operator from the LHC experiment. For this purpose, we implement 
the interactions described in section [2] in Madgraph5 with the aid of Feynrules [21] 
and perform MC simulations for three channels as discussed above. 

In order to put constraints on the operators involving light quarks (up, down, charm 
and strange quarks), we perform the mono jet analysis. We find the optimized 
cut on $ T by maximizing the ratio of the efficiency for accepting the signal to the 
upper bound on the cross section times the acceptance for the non-standard model 
events. For the reference point to obtain the optimized value, we choose the dark 
matter mass of 100 GeV and the interaction with up quarkj^ The results are sum- 
marized in table HJ We find that stronger cut on $ T gives more severe bound on 
the interactions for the Majorana fermion dark matter. It is possible that cut on $ T 
larger than the value we have shown in sect ion |4 . 1 . 1 1 gives more severe bound. How- 
ever, stronger cut will lead to larger systematic error thus we do not consider this 
possible improvement in this paper. The current and future expected constraints on 

lr rhc optimized value do not change more than 50 GeV even if we choose other reference points. 
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Real scalar 


Major ana fermion 


7 TeV 


350 GeV 


400 GeV 


8 TeV 


500 GeV 


500 GeV 



Table 4: Optimized cut on $t for the real scalar and the Majorana fermion dark matters. 

each operator are shown in the first and the second columns of figure E] and H] in the 
next section. 

In order to put constraints on the operators involving bottom and top quarks, 
we perform the mono b-jet +$t an d the top quarks(s) +$? analysis, respectively. 
For the top quarks(s) +$t analysis, the optimized cut on Mft turns out to be 900 
GeV for both the real scalar and the Majorana fermion dark matters. The future 
expected constraints on the operators are shown in the third columns of figures |3j 
andU 

5 Results 

In this section, we discuss limits on the dark matter effective operators mentioned 
in section |2j We have already explained in the previous two sections that the limits 
come mainly from WMAP, Fermi-LAT, XENON100 and LHC experiments. Since the 
structures of the operators are very different for the real scalar and the Majorana 
fermionic dark matters, we discuss the results of these two classes of interactions 
separately in following subsections. In all cases, we vary dark matter mass from 
10 GeV to 1 TeV and calculate the limits and future prospects in terms of the 
combination of the cut-off scale (A) and the coupling (c), yfcj A. 

5.1 Scalar dark matter 

Let us first discuss the (expected) limits for the real scalar dark matter interactions. 
All results are summarized in figure |3J where limits from cosmology, direct detection, 
and indirect detection experiments are shown as brown, green, and blue lines, while 
the current limit and future prospect from the LHC experiment are shown as magenta 
and pink lines, respectively. It can be seen that the limits obtained from the indirect 
detection are much stronger than those from the LHC experiment except for the 
operators <f) 2 it and (jP't^H. This is because the production cross section at the LHC 
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experiment is small for the scalar dark matter, while its annihilation cross section is 
large as there is neither helicity nor p-wave suppression. 

The limit from the direct detection experiment is also very strong for the scalar 
interactions, since they contribute to the spin-independent scattering cross section. 
The limit on the operator 2 cc is stronger than that on cf) 2 ss. The reason is the 
following. The hadron matrix element (iV|ss|./V) is small (/y s ~ 0.009), while the 
interaction with c quark is expressed by the interaction with gluons at low energy 
and the matrix element (N\G* v G a » v \N) is large (f Ta ~ 0.9431). 

It is worth noting that the lines obtained from the cosmological limit are almost 
horizontal except the top quark mass threshold for the operators <j) 2 it and (pH^t. 
This is because the annihilation cross section of the scalar dark matter is controlled 
by the dimension-5 operator after the electroweak symmetry breaking, and the cross 
section depends very weakly on the dark matter mass. 

It is also clear that the scalar interactions with u and d quarks are not favored 
from current data, whereas the interaction with top quark is the most difficult one 
to constrain. Even with 100 fb _1 data at 14 TeV run, it would be challenging to 
cover the most of the parameter space consistent with WMAP data. 

5.2 Fermion dark matter 

We next discuss the (expected) limits for the Majorana fermion dark matter. All 
results are summarized in figure HI where the limits are shown in the same way as 
the scalar dark matter case. A remarkable point is that the limit from the indirect 
detection experiment is weak due to the helicity suppression. It is to be noted that 
the limit from the direct detection experiment is also very weak (and not shown in the 
figure), because the operators do not contribute to the spin-independent scattering 
cross section. It is therefore important to investigate the experimental way to search 
for the Majorana fermion dark matter at the LHC experiment. 

The mono jet +$t search gives strong limits on the dark matter interactions 
with light quarks (up, down, charm and strange quarks) in comparison with the 
limits obtained from the indirect detection experiment in the most of the parameter 
space, and so will the mono b-jet + $ T search on the interaction with bottom quark. 
Unfortunately, the limit from the 8 TeV run will not be much stronger than the one 
obtained in the 7 TeV run, since the number of SM background events increases 
as well as that of the signal events does. On the other hand, the top quark(s) + 
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$ T search will give comparable limits to those from the indirect detection for the 
interaction with top quark, which is a sharp contrast to the case of the real scalar 
dark matter. 

6 Summary 

We systematically investigate dark matter interactions with various quarks in the 
framework of effective operators. We consider the dark matter to be a singlet scalar 
or a Majorana fermion and thoroughly study bounds and future prospects of the 
dark matter discoveries from different astrophysical, cosmological and collider ex- 
periments. We have shown that the mono b-jet + p T , and the top quark(s) + $ T 
channels along with the mono jet + {$t) channel could be very important to search 
for the dark matter interacting with quarks. Here, we derive model-independent up- 
per limits on the cross section times acceptance for these three channels mentioned 
above at 7 TeV (8 or 14 TeV) run of the LHC experiment. It is to be noted that the 
model independent bounds derived in our paper can also be used to constrain other 
dark matter or new physics models with similar type of signatures. 

For the scalar dark matter with spin-independent interactions, the direct detec- 
tion experiment supersedes the LHC bound because of smallness of the production 
cross section at the LHC experiment. In all cases whether it is spin- independent 
interaction or not, the indirect detection limit from Fermi-LAT data on the scalar 
dark matter is found to be stronger than the LHC bound. However, in case of the 
Majorana dark matter interactions considered in our work, the LHC experiment 
plays very important role for detection or exclusion of such possibility as the spin- 
independent cross section is zero. Possibility of detecting the top quark interaction 
with the dark matter is unfortunately not so promising at the LHC experiment even 
with 100 fb" 1 data at 14 TeV run. 
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Figure 3: Limits on the operators describing the interactions between the real scalar 
dark matter (f> and various quarks. Those from current LHC 'LHC(7TeV)', future LHC 
'LHC(8TeV) & LHC (14TeV)', cosmological 'Cosmology', direct detection 'Direct', and 
indirect detection 'Indirect' experiments are shown as magenta, pink, brown, green, and 
blue lines, respectively in each panel of the figure. 
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Figure 4: Limits on the operators describing the interaction between the Majorana fermion 
dark matter x an d various quarks. Those from current LHC 'LHC(7TeV)', future LHC 
'LHC(8TeV) k LHC (14TeV)', cosmological 'Cosmology', direct detection 'Direct', and 
indirect detection 'Indirect' experiments are shown as magenta, pink, brown, green, and 
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